(NASA-Ttt-83949) AN APPLICATION OP TBE 
BfiAONBfiCK flfiTflOD TO TBE MAaGI*RUuINCWlC2 
FIELD REPRESENTATION (NASA) 24 p 
HC A02/MP AOl CSCL 


G3/74 


Ne2-31108 

Unclds 

30365 


IVIASA 

Technical Memorandum 83949 


AN APPLICATION OF THE 
BRAUNBECK METHOD TO 
THE MAGGI-RUBINOWICZ 
FIELD REPRESENTATION 


Robert Meneghini 


JUNE 1982 


National Aeronautics and 
Space Administration 

Goddard Space Flight Center 

Greenbelt, Maryland 20771 



AN APPLICATION OF THE BRAUNBEK METHOD TO THE 
MAGGI-RUBINOWICZ FIELD REPRESENTATION 


Robert Meneghini 


June l')82 


GODDARD SPACE FLIGHT CENTER 
Grcenhelf. Miir> laiul 




AN APPLICATION OF THE BRAUNBEK METHOD TO THE 
MAGGI-RUBINOWICZ FIELD REPRESENTATION 





ABSTRACT 

The Braunbek method is applied to the generalized vector potential associated with the Maggi- 
Rubinowicz representation. Under certain approximations, an asymptotic evaluation of the vector 
potential is obtained. For observation points away from caustics or shadow boundaries, the field 
derived from this quantity is the same as that determined from the geometrical theory of diffraction 
on a singly diffracted edge ray. ITie paper concludes with an evaluation of the field for the simple 
case of a plane wave normally incident on a circular aperture, showing that the field predicted by 
the Maggi-Rubinowicz theory is continuous across the shadow boundary. 
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AN APPLICATION OF THE BRAUNBEK METHOD TO THE 
MAGGI-RUBINOWICZ FIELD REPRESENTATION 

INTRODUCTION 

The Maggi-Rubinowicz (M-R) representational’^) for the field is derived by means of Stokes’ 
theorem which transforms an open surface integral, the integrand of which has zero divergence, mto 
a line integral around the boundary of the surface and certain singular contributions interior to this 
boundary. This type of field representation is usually appUed under a physical optics approxima- 
tion, which, in the case of diffraction by an aperture, is taken to denote the approximation by 
which the actual field within the aperture is replaced by the unperturbed incident field. There is 
nothing in the theory, however, that requires this assumption. 

Another common high frequency approximation is the Braunbek method^^'*’^), which takes the field 
near an edge point as that which would exist if the edge were replaced by a half plane oriented in 
the plane of the aperture with its edge lying along the tangent vector at the point of interest. An 
application of this method to the various surface integral representations for the field has been 
shown to produce results that closely resemble those predicted by the geometrical theory of diffrac- 
tion (GTD)^^'^). Since the M-R theory is, in a sense, only a restatement of certain surface integral 
representations, there is reason to expect that the Braunbek method should also serve to illuminate 
the relationships between the M-R theory and GTP 

In this paper, the Braunbek method is appUed to the generalized vector potential of M-R^'). Under 
certain approximations, a closed form expression for this quantity is obtained. With this vector 
potential, an asymptotic evaluation of the field is carried out for observation points away from 
caustics and shadow boundaries. It is found that the expression for the field is identical to a result 
derived by Keller et al.^^). In the last section of the paper, an explicit calculation of the M-R field 
is carried out for a simple case, showing that the field is continuous across the shadow boundary. 
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THE RAYLEIGH INTEGRAL REPRESENTATION 

Consider a thin screen in the z = 0 plane with an aperture A. The tangent to the rim of the aperture 
is taken to be continuous. If the sources of the incident field are located in the z < 0 half space, 
and if the field is zero on the screen, then the Rayleigh representation for the field in the z > 0 half 
space is^®^ 

u(x)“2ju(x')v'G* zds' (1) 

A 

The region of integration is the aperture; x,x' are, respectively, the position vectors of the observer 
and of a "source” point in the aperture. The primed derivatives are to be taken with respect to the 
“source coordinates”. The quantity G is the Green’s function, where, for an e“‘‘’^* time convention 

G “ (2) 

and 

V'G»-f(ik-l/r)G (3) 

f »(x-x')/| x-x' I (4) 

r»|x-x'| (5) 

The quantity k in (2) and (3) is the wavenumber which equals lufk where X is the wavelength of the 
incident field. 

The M-R analogue of ( 1) is^®^ 

u(x)*#(W' W|) • 8d£' + ^ • 8d£' (6) 

i 

The path of integration of the first term, the boundary diffraction integral, is taken in a counter- 
clockwise sense (when viewed from z > 0) around the rim of the aperture. The second term repre- 
sents the integrations around the singularities of W, taken in a clockwise sense. Following Miya- 
moto and Wolt^ ' ^ the singularities of W are assumed to be isolated poles, finite in number. A 
corresponding term that accounts for the singularities in Wj is absent in (6) because W| is never 
singular over the aperture^ 
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The generalized vector potentials for the observation point and its image, W, Wj, are given 

m 

W = - G s X J e*^* V' u(x' + ts) dt + (7) 

0 

•• 

Wj * - G Sj X e^^^ V' u(x' + 1 Sj)dt + (W^)j (8) 

o 

where 


A A 

s = -r 

S| = - fj ■ s - 2(s • z) z 


(9) 

( 10 ) 


The end point contributions to the vector potentials, and (W^)j are defined by Miyamoto and 
Wolf^'\ In this paper they will be neglected. The effect these additional terms might have on the 
solution is discussed later on. 

SL'ce the Braunbek method is valid only in the vicinity of the edge, it is applied only to the first 
term of (6). For the second term of (6), the vector potential corresponding to the physical optics 
approximation will be used. 

APPLICATION OF THE BRAUNBEK METHOD TO THE VECTOR POTENTIAL 
For aperture diffraction, the physical optics approximation takes the field in the aperture to be 
equal to the unperturbed incident field. For the same problem, the Braunbek method assumes the 
field near a point 2 on the aperture rim as that which would exist by replacing the edge by a half 
plane oriented in the plane of the screen with its edge aligned along the tangent vector at 2. 


In this paper, the incident field, u^, is approximated in the vicinity of an edge point 2 by 
Uq = A(2) The phase \//(2) is assumed to satisfy the eikonal equation 

VV*V>»1 (11) 

so that a unit vector p can be defined by v V ® P- Again, following Keller et al.^^\ the incident 
field at 2 is associated with a plane wave propagating the direction p with an amplitude A{2). 


The solution to this half plane problem, with u = 0. on the plane, can be written in the form 
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(Appendix) 


u = u* “ u' 


( 12 ) 


u. 

01 




. I 

y‘r a sin®o'- z COS0 q 1 p («.) 

I 2 

«1 =-e‘7i 
02--€'72 



i 


= sgn (cos H(0?0o)) 


sgn (x) = 


+1 x>0 
-1 x<0 


F 



X 



dr 


The notation follows that of James^^^ and is defined in the Appendix. 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


To use ( 1 2) through ( 1 7) in (7) and (8) it is convenient to first express the above formulas in terms 
of quantities independent of a particular coordinate system. At each point 8 on the rim, we choose 
the orthogonal vectors (8, n, z) where 2 is the tangential unit vector of the rim at 8 (in a counter- 
clockwise sense when viewed from z > 0) and n is the unit vector from 8 directed into the aperture. 
Defining | to be a unit vector from 8 to an arbitrary point and t the associated distance along the 
vector to this point, then 


1 . 2 = - sin6 sin0 

(19) 

1 • h = - sin6 cos0 

(20) 

1 • 6 * cos6 

(21) 

p • z * sin0Q sin0Q 

(22) 

p • n = sin®Q COS0Q 

(23) 

p • 8 = cos6q 

(24) 

p = t sinfl 

(25) 

z = t cos6 

(26) 
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where 

P = V' ^ 

(27) 

Expressing the variables in (12) through (17) in terms of the quantities p, i ft, z, t, then 



ui(|, ,) = ir-» A < P 'i' t F, ««i) 

(28) 



(29) 

where 

1= ee + nft + zz 

(30) 


1-= 58 + fin - zz 

(31) 


(l.t) = e‘ (l)7i (it) 

(32) 


tt2(it)=«^(i 72 (it) 

(33) 


7l (i t) = (kt)^ (T(l)- P - It • 

(34) 


72 (it) = (kt)^ (T(t)-p-lt-i^ 

(35) 


£' (1) = (sgn (i • z) «i -i, <» • 1. • ♦ P • 

(36) 


±(({•1, • l)“+i • 


with 

• 




(37) 



(38) 


i “i- 

(39) 


To compute the V 'u appearing in the integral of (7) and (8). a spherical coordinate system is chosen 
with the origin at a point i on the rim with the polar axis in the ? direction. In this sytem. the a g 
ments of v'u. x' + 1 Und x' + t s, can be written as t| with { evaluated at s and s,. respectively. 
However, we will continue to write the arguments of u as t since these quantities do not 
always appear together. In spherical coordinates 


(401 


I' ' t 


PI)N, . ' 
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where 6 is the polar angle measured from the 8 vector to k and 0 is the azimuthal angle measured 
between the -n direction and the projection of | onto the n-z plane. For an observer on the 8 
axis, a positive <f> angle is measured in a counterclockwise sense from the -n axis.; Note that these 
are the same angles as defined by ( 19H21). Equation (6) shows that only the 8 component of 
(W “ Wj) is needed to compute the field. Since |x|*0, |xd *8®0 and | X 0 • i * sin0, then 
for the vector potential W , 


so that 


, ^ . a I 9«(i t) 

Itxv'ud.ol-t-Y-j; 


, ^ pe“^' 3u(it), 

V/-8-GJ-— 

A similar simplification can be carried out for W|. Therefore, 


(41) 


(42) 



A gikt / 3 „ . 3u £ V 1 \ 

(43) 

a 

In computing — u 

, the following relationships are used 

(44) 



(45) 



(46) 



(47) 

which give 



30 " 

A - ^/4) . j ikt [ B(|) («J ) -C(|) F,(a.) 

1 (48) 


.Ki 


Mill ekt m + iffi gi(|) ~ ^ e''“ *“2 e^({) 


y\ 




with 


P(5) = P -i* 5 


(49) 
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i^oo\ ^rr 


Q(l) = p • i • i 
c(|) = p • (L X fi) • I 


In finding the integrand of (43) from (48), a considerable amount of simplification is possible by 
using the following relationships between s and Sj 

p‘i*s*p‘i*§|;p*l*lj = p‘i*s (f 

p • (Jj X g) • s = -p ‘ (It X g) • spp • X g) • §i = -p • (|j X g) • s a 

eW§) = -e'(§i) (f 

e‘(si) = -e'(s) (S 

«i (s,t) = -0(2 (sj.t) (5 

«j (si,t) = -a2 (s,t) (5 


Then (43) becomes 


(W-Wi)*g* 


0 M ^ A 

“ C(s) e^^ F(a 2 )l dt + k J e’*^^ e^* e‘(s) 

o L 

-^^2ctQ(sO + ia^gf(5)jdt| 


, ^ Q(s) 


where 


F(x) » C e‘^^ dt (60) 

P(s), Q(s) B(s), C(s), are given by (49)-(53) with S replacing In (59) the arguments of «| and «2 
are s and t. Therefore «j and 02 can be evaluated from (32)-(36) by replacing $ with I. 

It is worth mentioning that the application of the Braunbek method here is somewhat different 
than in Keller et al.^®^ where the method was used to obtain an approximation of the field near the 
edge in the plane of the aperture. Here, the method is also used to approximate the field near the 
edge, but the vector potential requires that the field be evaluated along the path of the vectors s 
and S|. These vectors will be in the z = 0 plane only tor observation points on the aperture or 


• 'i m nr-ff ny fsi 
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screen. One consequence of this difference is the retention of the and quantities in (59). In 
Keller et = 1 and = -1, irrespective of the location of the observer. 

To find a closed form expression for the vector potential (W - W|) • fi, (59) is evaluated by using 
the first term in an asymptotic expansion for large k. Since the Braunbek approximation is valid 
only in the vicinity of an edge (i.e. near t » 0), the upper limit in the integrals of (59) will be 
ignored. In fact, ignoring the upper limit is more than just a matter of convenience since the first 
two terms of (59) do not converge. This problem is not necessarily a deficiency of the M-R theory 
but perhaps a consequence of having ignored the and (W»)i contributions to vector-potential. 
The results obtained by Miyamoto and Wolf^*^ for a convergent spherical wave (which is similar to 
the Braunbek method in the sense that neither satisfies the Sommerfeld radiation condition in 
three dimensions) suggest the possibility that the full expression for the vector potential would 
converge and even admit an evaluation of higher order terms in the asymptotic series. Since we 
have not been able to compute W. and (W«)j, we obtain an expression for the first two terms by an 
integration by parts, evaluating the result only at the lower limit, and using the formula^^* 


i” «'>"«' (kh'(0))» 


( 61 ) 


The last two integrals of (59) are of the same form as (6 1) and can be evaluated immediately. 

From the results of the calculations and simplifications, the vector potential can be written 

(W-W,)*£=^ (l+T(§) + (&-«)(s-fi))^ 

“ ‘ 4jrr 


e* (s) p ♦ ^ X g) ♦ § 

I (T(§)-0 •§)^ (l+^*i‘S)“ (T (S) 


(s)p-(;Lxe)-s 
-P-Xt *8)^(1 +&*i* S). 


(62) 


This formula can now be compared with the first tenn in the asymptotic expansion of (W - Wj) • £ 
under the physical optics approximation,^* 


(W-W,)-gJpo.= 


A(g) 


^ f sx p ♦ S _ s, X p • 5 ‘ 

J +8 • ^ 1 + Sj • ^ . 


47rr 


(63) 






onic?rv''^'' r 

OF PO ! t 


The comparison is made easier by noting that 

P‘^^xS)*§ = (§Xf>)‘S 
p • X g) • s = (§i X p) • S 
(1 +p * 1 * s) = 1 + P * s 


(i+p*i*s) = i+ p*sj 


where 


si = s- 2 (§*i)i OO) 

Therefore, the first terms in (62) and (63) differ by the factor (§) { 1 1 + T(S) + (&•«(§•«)]/ 

[T (s) -p • • §]|^; the second terms by the factored (s) |[1 +T (§) +(^ • 5) (§ • 5)]/[T (s) - 

p-i-f. 

COMPARISON OF M-R AND GTD 

As a check on the vector potential given by (62), the boundary diffraction term of ( 6 ), Ug, is evalu- 
ated for large k when the observer is not close to a shadow boundary or caustic, i.e. 

Ug(x) = #(W-W|)*«dK' ( 68 ) 

where (W - Wj) • 8 is given by (62). 

Recognizing that the saddle points of ( 68 ) occur when 

~( 0 +r) = O (69) 

then since 3 i/// 9 £' = 5 • and 3r/92‘ = s • £ the saddle points are given by points on the rim £ = £54 
for which § • p =» - S • s. Using the formula for an isolated first order saddle point gives( ^ ^ > 


Ub.(2^ 


(27t)^ (W - Wt) ♦ £ 


(kir + r"i)^ 


where 


+ ^ (71) 

* 9£'2 9£'2 

To compare (70) with the result of Keller et al.(^\ (W - W) • 8 is evaluated at 8 • p = - £ • s and the 
result is expressed in terms of the angles 0^, 0. 0^,, (* defined by ( 19) • (27). In making this trans- 
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rotation, i. i. imporun. .o note ,ha. ,ha vector | in ( 19) • (2 .) U talc™ to be f , the unit vector 
from the rim to the observer, where f - s. 


This procedure gives 


_ w \ . g — — -[e‘(s) sgn(siny2(0o“^^^ 

4ffr 


( 72 ) 


(W-Wi 

+ e^(s) sgn(sin V 2 (<I>q+<I>)) sec Vaf^-hAo)! 
where both sides are evaluated at 2 • p = “ 2 • s. 

Since only observation points in the z > 0 half space are considered, it is not difficult to show that 

e*(S) sgn {smVi(<t>Q~<t>)) “ “ 1 
e^(s) sgn (sin Vi(^g+^)) * 1 

a- ..e«r^ Kv if<.iier ftt al let (i„ ® a + ff/2 and 0 = 6 + it/2. 

To account for a different convention used by Keller et ai. , Vq 

Substituting these equations into (72) and using (73), (74) gives 

l AcW^.)e-/4tP(*‘^0 ^ 1 (75) 

Ub(x)“-\ 2(27rk.|<//"+r"l)^ 

This ia the aante expression obtained by appiyinj the Braunbek approximation directly to the Ray^ 
,ei,h representation and is also the field predicted by GTD on rays singly diffracted from an edge . 

asymptotic evaluation of the m-r representation for a simple case 

Although we have been unable to show in general that the vector potential given by (62) leads to 
fields that are continuous across the shadow boundary, there are special cases where this continuity 
can be shown explicitly. It should be mentioned that under a physical optics approximation, Otis, 
et al.( 12) and Takenaka et al.( '^l have obtained asymptotic evaluations of the M-R representation 
for an on-mtis gaussian beam incident on a circular aperture. Their work has shown the field to be 
everywhere well behaved. Here we will treat the simpler case of a plane wave nonnally Incident on 
a circular aperture of radius a, but use the vector potential given by (6’), 

For a normally incident plane wave of unit amplitude. (62) reduces to 
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ikr 


(W~W.) • = +(f •!, • • 

* Am * 


(f X fe * S) (-r) 


((f * * r)^ + r • i)*^ ( 1 - f • 2) 


(f, X z • e) (-f) 


((f • i, • f)^ + f[ • z)'‘^ ( 1 - fj • 2) 

where s = -f and f j = f - 2 (f • z) z. 


(76) 


The total field in the x-z plane is 

u(x) » # (W - Wj) • « dfi' + [U(x + a) - U(x - a)] (77) 

where the second term represents the geometrical optics field. The boundary diffraction integral 
can be written as the sum of I j and I2 , 

Ij +l2«^(W-Wj)*fid£' (78) 


where 


and 


3ff/2 

(79) 

3ff/2 * V * 1 

(80) 

II 

(81) 

^ €*(-r)(l +(rMj • r)^)^ 

/l(0 1 “ ^ 

r ((f • * f)” + f • z)” 

(82) 

AW) 

r((f *J[j • f)^ +fj • z)^ 

(83) 

eVf) = 1 

(84) 


The expression for It can be evaluated by means of the standard formula for an isolated saddle 
point. The saddle points are determined from the equation C • ^'r = - $ • r = 0. Since ? = - x sin0' 
+ y cos(f' and f = (x(x-a cos0') - y a sin0' + zz)/r then saddle points occur at 0. =0 and 0, = it. 
Noting that e^[- r (0jj^ )| = sgn (a-x) and H0j^)l = sgn (a + x). a straightforward calculation 
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gives 

^ i /l.'\ Q^rittff/4 [a-hxl 
2 \ity |(x (tj-z))^ (fj+z) 

with 

fj = ((a - x)^ + 

T 2 “ ((a + x)2 + 

By expanding (85) about the small parameter xTa, it can be shown that I 2 is continuous across the 
shadow boundaries at x * ± a. 

Turning to the expression for Ij , it can be shwn that the steepest descent path may be close to the 
pole at r • 2 = 1. Therefore, the standard formula no longer applies^* In the following develop- 
ment we make use of the fact that the field is symmetric about the z-axis. Therefore, it is sufficient 
to consider only observation points in the x-z plane for which x > 0. Since the saddle points of 
are th«! same as for I 2 , then 0^^ * 0 and = ir. In the complex 0-plane, the steepest descent path 
(SDP) which passes through 0jj is given by Im (i r (0)) * Im (i r (0jj)) where Im denotes the 
imaginary part. Graphs of the steepest descent peths are given by Takenaka et al* 


gikr2“lJr/4 u- 


(x (r 2 “z)) 


|a-x| 

(r2+z) 


(85) 


z2)H 

z2)H 


(86) 

(87) 


Solving the equation r • 2 * 1 for 0 yields poles at 


0p^ » i Cosh“^ q 

(88) 

0_ » - i Cosh“* q 

(89) 

where 

q = (x^ + a^)/2ax 

(90) 

Because of the presence of r in the denominator of the integral, I j 

also has algebraic singularities at 

0^j » ± i Cosh q' 

(91) 

q' = (X' + a^ + z“)/2ax 

("2) 


For z » a, these contributions can be ignored. 
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In the deformation from the original integration patii to the steepest descent paths, no poles are 
crossed so the contributions from the residues that oceur for the case of a gaussian beam^**^^ are 
absent in this case. Furthermore, the poles can be located only in the vicinity of the SDPj ( since 
we have chosen x > 0). The integration along SOP, therefore can be evaluated in the same manner 
as l2- 


Using the formula derived by Felsen and Marcuvitz^ * ^ ^ (with a modification to account for the 
presence of two poles) then 

Ij ^ “ a |^sgn(lm bj)2lR57T^^ e”*^ Q(~i bj sgn (Ini bj)) + * 


where 



(93) 


Q(x) e"‘^ dt (94) 

X 

The quantity Rj is the residue of /| (0)(fX2*^)/( 1-f ‘z) at the pole 0pj (i = 1 , 2) and 

hj = (-2/i|r"(0^)!)^ i=l,2 (95) 

t'i=lUr(0,|)-r(0pj))l^ i=l,2 (96) 

where the argument of bj is to be chosen the same as the argument of the quantity (0pj - 0^^ )/h| 
as 0pj approaches 0j^ ^ * \ 


To compute Rj. the following result is used^***^ for the residue of t(0)/p(0) 



as X approaches a, i.e„ as the observer approaches the shadow boundary, 1,^ + 1 1 ^ = 0 and as the 
observer crosses the shadow boundary, the discontinuity in exactly cancels the discontinuity in 
the geometrical optics field. Since I 2 and the remaining term of Ij are continuous, this implies that 
the total field is continuous across the shadow boundary. For points away from the shadow 
boundary (k^ tj » 1). + Ijb is o(k“^) and the I 2 and the remaining terms of yield the 

standard asymptotic result. 

Points on the z-axis, i.e. the axis passing through the center of the aperture, correspond to caustics 
of the diffracted field and equation (106) is not valid on or near such points. For such points, the 
procedure described by Takenaka et al.( ‘3) could be followed. However, we will not pursue the 
problem here. 




CONCLUSIONS 

In an attempt to establish a correspondence between the M-R and GTD theories, an approximate 
form of the vector potential has been derived through the use of the Braunbek method. For the 
aperture diffraction problem, the two theories predict similar results for observation points away 
from caustics and shadow boundaries. To examine the behavior cf the M-R results near the shadow 
boundary an explicit calculation was performed for the simple case of a plane wave incident upon a 
circular aperture. That the M-R field was continuous across the shadow boundary in this case, 
suggests that the vector potential given by ( 62 ) might yield reliable results near a shadow boundary 
for a more general incident wave. 
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Appendix: Diffraction from a half-plane 


If an incident plane wave of the form 

ujj = A cos(i^^ ) sinflo g-jkz cosBq j 

is incident on a half plane the edge of which is aligned with the z-axis and oriented in the x-z plane 
(x > 0), and if the total field, u, on the screen is zero, then the field at a point x is 

u = u‘ - (A2) 

where 

u‘ = U(e‘) ujj + uj, (A31 

uf = U(ef) (A4) 

and 


“d ® 


I I 

- A K_(t:kpsin0^)^ | cos-(«4«^,) |) + 2 cost?o) 

o 

' ■ 1 

6^ = sgn [ cos-(0T0q)I 

A* 

K_ (y) = I e“i‘^ dt 


.. ( • x> 

U(x)= ^ ^ 

( 0 X < 


(A5) 

(A6) 

(A7) 

(AS) 


An el‘^' time convention is used to obtain these results. Apart from some sliglit modifications - 
in particular, the addition ot the factor A which provides an amplitude and phase reference 
for the Braunbek method - the above formulas were obtained from James^‘*\ 


From equation (A 1 ), the direction of the incident wave, p, is 

p = “ cos0^^, sin6^ X - sin0^, sind^, y + i (A'J) 

The a/imuthal angle 0^, is measured in a counterclockwise sense from the x-axis to the projection of 
the vector ~p onto the x-y plane. The polar angle is measured between the z and p vectors. On 
the other h.ind. the coordinates of .m observation point x with respect to an origin on the edge of 
the half-plane is 


X = p COS0 X + p sin0 y + /> 


(AID) 




fij. '^»; •’' 


where 


p = I X I sin0 
z = I X I cos6 


(All) 
(A 12) 


The angle 0 is measured in a counterclockwise sense from the x-axis to the projection of the vector 
X onto the x-y plane; 6 is the angle between z and x. 


From (A3) and (A7) 

u‘ = A e'^'^ co^<H>o) “ z c°s®o] 

where is given by the argument of K_ in (A5) usi 


U(e‘) - 7 t“^ e‘ dt 

>1 J 


(A 13) 


ng the upper sign. 


Since 


m 

1 = y e"7^^ dt 

the expression in the brackets of (A 13) can be rewritten as 

M 

tT^ y dt 


where 


so that 


«l=-e‘7i 


Similarly 


m 

a ^ g-jk^ + jff/4 gjk(pco<^o) an ®0 “ 2 <»s®o] ^ dt 


•• 

yf _ ^ + jrr/4 gjk[pcos(^o) an^o “ z cotfo) ^ dt 


where 


a-, * - 7-» 


and where 7 -. is the argument of K_ in (A5) using the lower sign. 


(A 14) 


(A 15) 


(A16) 


(A 17) 


(A 18) 


(A 19) 
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To account for the convention used in the text, we substitute j = -i into (A 1 7) and (A 1 8). 
This gives equations (13) of the text. 
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